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1. Introduction 

Security of sustainable energy is one of the most important con¬ 
cerns for world societies [1]. A number of commercial enterprises 
are required for supplying climate-friendly energy technologies 
in future [2-5], including electricity produced from solar and wind 
energy [6-9 , bio-fuels, and carbon dioxide (C0 2 ) capture and stor¬ 
age (CCS) [10-15 . In this regard, scrutinized planning with respect 
to the prospect plays an essential role toward better future. Green¬ 
house gases emission is one of the restrictions of utilizing fossil 
fuels which encourages development of higher efficiency and max¬ 
imal electric [13-16[ in energy production systems. Use of Stirling 
engines results in green house gases emissions reduction into the 
atmosphere. The Stirling engine has a high potential to be applied 
for converting heat into the mechanical work with a high thermal 
efficiency. In theory its thermal efficiency might be as high as the 
Carnot efficiency when it is ideal regenerator, reversible, and hav¬ 
ing isothermal compression and expansion processes. Stirling en¬ 
gine is an external combustion engine, so it can be powered by 
various heat sources and waste heat [17-24 . Blank et al. [20] stud¬ 
ied the power optimization of an endoreversible Stirling cycle and 
provided an estimate of potential performance for a real engine. 
Thombare and Verma [21] gathered the available technologies 
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and obtained achievements with regard to the analysis of Stirling 
engines and, at the end, presented some suggestions for their 
applications. Formosa and Despesse [23] conducted the modeling 
by means of the isotherm model in order to investigate the effects 
of dead volumes on the engine’s output power and efficiency. 

The Stirling engine cycle can be defined as a closed regenerative 
thermodynamic cycle consists of compression and expansion of the 
working fluid at different temperatures [25—35 . Tlili et al. [33] devel¬ 
oped a first law model with additional consideration for an internal 
irreversibility parameter to examine the effects of irreversibility on 
net work, heat addition, and thermal efficiency. Their irreversibility 
parameter originates from the second law of thermodynamics for a 
real cycle. Their results regarding dead volume and regeneration ef¬ 
fects support the observations of Kongtragool and Wongwises [35]. 
They also reported that larger irreversibilities in the Stirling cycle 
thermodynamic processes lead to smaller maximum network output 
and thermal efficiency. However, they did not quantify the roles of 
the heat exchangers, regenerator, and dead volume irreversibilities 
on the overall system irreversibility. 

The enhancement of finite-time thermodynamics [36-40 , lays 
on a useful tool for performance analysis of real engineering prob¬ 
lems. The finite-time thermodynamic performance of the Stirling 
engine has been scrutinized by a great number of investigators 
[39-46]. The effect of heat transfer on the power output of a Stirling 
engine has been inquired by many researchers 4,23-30]. An optimal 
efficiency for maximum power output of a solar Stirling engine con¬ 
sidering imperfect regeneration is achieved by Petrescu et al. 
[39,41]. Moreover, finite-time thermodynamics analysis of heat 
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Nomenclature 



Ah, A l 

heat transfer surface area between the endoreversible 

t 2 

temperature of the working fluid during the isothermal 


Stirling engine and either the high-temperature source 


external heat rejection process 


or the low-temperature sink 

Th. Tl 

temperature of the heat source and the heat sink, 

Areg 

heat transfer surface area of the Stirling engine regener- 


respectively 


ator 

t c 

time required to accomplish the isothermal compres- 

Cv 

constant volume specific heat 


sion 

Cl, 

decision index in TOPSIS method 

th 

time required for the heat transfer process from the 

d 

deviation index 


regenerator to the working fluid 

d i+ 

deviation or distance of ith solution from ideal solution 

tR 

time spent on the regenerative processes 

di- 

deviation or distance of ith solution from non-ideal 

i 

ith objective 


solution 

j 

jth solution 

m 

mass of the gas 

U H .U L 

Overall heat transfer coefficient between the engine and 

P 

power output 


either the low-temperature sink or the high-tempera- 

0 . 1-2 

isothermal heat transfer from the engine to the low- 


ture source, respectively 


temperature sink 

U re g 

overall heat transfer coefficient of the regenerator 

02-3 

heat transfer from the regenerator to the working fluid 

X 

vector of decision variables 

03-4 

isothermal heat transfer into the engine from the high- 




temperature source 

Greek letter 

0 . 4-1 

heat transfer captured by the regenerator from the 

A 

ratio of volume during the regenerative processes (vol- 


working fluid 


umetric ratio) 

R 

gas constant 

n 

efficiency 

T, 

temperature of the working fluid during the isothermal 

ILL 

fuzzy membership function 


external heat input process 

3 

the Stefan’s constant 


engines is often limited to systems having linear heat transfer law 
which depends on the temperature differential both the reservoirs 
and engine working fluids [22,44-50 . Yagi et al. developed a math¬ 
ematical model for the overall thermal efficiency of solar powered 
high temperature differential dish Stirling engine with finite heat 
transfer and irreversibility of regenerator and optimized the absor¬ 
ber temperature and corresponding thermal efficiency 42]. Tlili 
investigated effects of regenerating effectiveness and heat capaci¬ 
tance rate of external fluids in heat source/sink at maximum power 
and efficiency [44]. Kaushik and Kumar studied the effects of irrever¬ 
sibilities of regeneration and heat transfer of heat/sink sources [46]. 

Solution of the multi-objective optimization problems is an ex¬ 
tremely difficult goal which requires the simultaneous satisfaction 
of a number of different and even conflicting objectives. Evolution¬ 
ary algorithms (EA) were initially extended and employed during 
the mid-eighties in an attempt to stochastically solve problems 
of this generic class [47 . A reasonable solution to a multi-objective 
problem is to investigate a set of solutions, each of which satisfies 
the objectives at an acceptable level without being dominated by 
any other solution 48]. Multi-objective optimization issues illus¬ 
trate a possibly uncountable set of solutions which is named Pareto 
frontier, whose evaluated vectors depict the best possible trade¬ 
offs in the objective function space. In turn, multi-objective optimi¬ 
zation of various thermodynamic and energy systems have been 
concentrated recently 49-56]. 

In real model on Stirling engines, heat regeneration is supposed 
be considered imperfect. The optimal power theory design of the 
Stirling engine may be complex in this case. Assuming an ideal 
Stirling engine can lead to acceptable results in power output 
and thermal efficiency. This article concentrates on these elements 
using two scenarios. 

2. Scenarios 

2.2. Thermodynamic modeling 

The regenerative endoreversible Stirling cycle is depicted in 
Figs. 1 and 2. Stirling cycle has two irreversible isothermal 


processes and two reversible isobar processes. External heat trans¬ 
fer processes are carried out in finite time. Process 1-2 is per¬ 
formed at the constant temperature T h and process 3-4 is 
performed at the constant temperature T c . 

In this work, it is assumed that a perfect regeneration is 
occurred in the regenerator. 

2.2.2. Mathematical model 

At a perfect regeneration, the amount of rejected heat from the 
hot working fluid to the regenerator mesh during processl-4 is 
equal to the amount of absorbed heat by the cold working fluid 



—► 

S 


Fig. 1 . Temperature-entropy diagram of the regenerative endoreversible Stirling 
engine. 
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Fig. 2. Functional schematic of Stirling engine (note that the opposing pistons are 
90 degrees out of phase). 

during process 2-3. This assumption is acceptable, as regenerators 
with effectiveness of 98% and 99% have been reported [17,20]. 


Therefore, we have; 


10-2-3 1 — 0-4-1 

a) 

Q-h = U h A h (Th - Th)th = mRThLn(A) 

(2) 

Q l = U l A l (T c - T L )t c = mRT c Ln(X) 

(3) 


Thus, it can be shown that the net work of the cycle is obtained as 
follow [17,20,42-46 ; 


W net = mR(T h - T c )Ln(A) 


(4) 


where m is the mass of the working fluid, R is the gas constant, T is 
the temperature of working fluid in each state of the cycle as de¬ 
scribed by subscripts and 2 is the compression ratio. The output 
power is obtained by dividing the net work to the time period of 
the cycle [17,20,42-46]; 


P = 


W net mR(T h - T c )Ln(2 ) 


t 


th + fc + t R 


(5) 


where t R is the total time spent for the regeneration processes and t h 
and t c are the times spent for isothermal expansion and compres¬ 
sion processes respectively. These time periods are defined as 
follows; 


th 


mRT h Ln(X) 

U h Ah(T h - T h ) 


( 6 ) 


where U H and A H are the overall heat transfer coefficient and the 
heat transfer surface area for the hot side heat exchanger, respec¬ 
tively. Similarly, U L , A L and U reg , A reg , similar values for cold side heat 
exchanger and regenerator. T H and T L are the temperature of heat 
source and heat sink respectively. C v is the specific heat at constant 
volume in kj kg -1 K -1 and is assumed to be constant. 

Assuming the ideal regenerative, ( LMTD) reg = 1 K. [4 . Substitu¬ 
tion of Eqs. (6)-(8) into Eq. (5) leads to the following expression: 


P= 


T h 


Tr 


2 C 


v 


_UhAh(Th— Th)(Tfi — T c ) UhAhi//Ar(T c — Tl)(Ti 2 — T c ) RU reg A reg Ln(X ) 


-i 


(9) 




Ar = 


U L 

U H 

Al 


A 


H 


For simplicity, we considered the x 
be written as follow; 


P = 


1 


and M 


x 


(9a) 

(9b) 

- Cv , so Eq. (9) can 


RLnA 


2 M 


i -l 


,^h^h(Xh — Xh)(l —x) U H A H i//A R (xTh-T l)( 1 —x) U reg A reg _ 


( 10 ) 


To determine the maximum power output, take the derivative of P 
with respect to Ti and obtain the optimum temperature of working 
fluid at state 1; 


dP 

So; 


= o 


VWrTl + xT h 


Ihopt x(CWr + C 

Therefore, the thermal efficiency is calculated as follow; 
P 

Substituting Eq. (11) into Eq. (10); 

1 x 


(ii) 


( 12 ) 


max 


2 M 


-i -i 


UhAh(Th— Tho P t)( 1 — x) U H A H i//A R (xTi lopt — Tl)( 1 — x) U reg A reg _ 


(13) 


For the maximum power output it should be considered U H A H = UA 
[4]. So Eq. (13) can be written as follow; 


P max — 


1 


X 


2 M 


-i -i 


(14) 


_UA(T H -T hopt )( 1 -x) UAi//A r (xT h op t - Tl)( 1 -x) U reg A reg _ 

Also, the corresponding maximum net work is obtained as follow; 

W ne L 0pt = mRT hopt ( 1 - x)Ln(X) (15) 

2A.2. First scenario results 

To have a numerical appreciation of the result, we consider the 
heat source temperature and cold source temperature 1300 I< and 
320 K respectively. Also U H A H = U L A L = 2 kW K“\ x = 0.5, 2 = 2, 
x/y = 1, A r = 1, and U reg A reg = 1000 kWK" 1 [20]. We consider stan¬ 
dard air as a working fluid with R = 0.287 kj K -1 kg~\ C v = 0.718 - 
kj K _1 kg 1 and m = 9.3 g, as a total mass of the working fluid. 


mRT c Ln(X) 

U l A l (T c -T l ) 

2 mc v (T h - T c ) 
UregAreg (LMTD) reg 


(7) 2A.2A. The effects of T H . As be shown in Fig. 3, in various UA , by 
increasing T H the output power increases. Also, in the certain T H , 
by increasing UA, output power increases too. As is clear in Fig. 4, 

(8) in various ( UA ) reg , by increasing T H the output power increases, also 
in certain T Ht with increasing ( UA ) reg , output power is increased. 
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The effects of temperature of heat source on the output power 
at various values of x are presented in Fig. 5. At temperatures above 
1100 K with increasing x, the power output is reduced. The slope of 
curve decreases by decreasing the difference between the temper¬ 
ature of working fluid in hot and cold spaces. According to Fig. 6, 
for different values of 2, with increasing T H the output power is 
increased. 

2.2.2.2. Effects of volumetric ratio (X). Also in Fig. 7, the effects of 
volumetric ratio on the output power at various UA are presented 
that have similar behavior to Fig. 6. 

Fig. 8 shows the effects of the volumetric ratio on the thermal 
efficiency at various values of temperature ratios, in which the effi¬ 
ciency increases by increasing 2 and as the temperature difference 
is decreased the thermal efficiency is reduced. Fig. 9 shows the ef¬ 
fects of volumetric ratio on the output power for different values of 
x that has similar behavior as the thermal efficiency be shown re¬ 
cently. Fig. 10 shows the effects of volumetric ratio on the output 
power for different values of ( UA) reg . 

2A23. Effect of UA. By increasing the UA of heat exchangers the 
output power increases rapidly in comparison with volumetric ra¬ 
tio and heat source temperature that means the heat exchanger 
and regenerator with high effectiveness coefficient are more effi¬ 
cient, as be shown in Fig. 11. 

2.2. Multi-objective optimization with evolutionary algorithms 

2.2A. General concepts of the multi-objective optimization 

Imagine a decision-maker that would like to optimize objec¬ 
tives which are non-commensurable. The decision-maker does 
not know the priority of the objectives. Moreover, all objectives 
are of the minimization type. For a minimization kind of objectives, 
it can be changed to the maximization type by multiplying nega¬ 
tive one. The procedure of a minimization multi-objective decision 
problem with objectives is represented as follows: 

Given a non-dimensional decision variable vector x = {xi,..., x„} 
in the solution space X find a vector x* that minimizes a given set of 
K objective functions/(x) = {fi(x*),... ,/&(**)}. The solution space X is 
generally restricted by a series of constraints, such as g/x*) = bj for 
j = 1,..., m, and bounds on the decision variables [48]. Overall, no 



Fig. 3. Variation of the power of the Stirling engine for different the heat source 
temperature and UA. 



Fig. 4. Variation of the power of the Stirling engine for different the heat source 
temperature and ( UA) reg . 



Fig. 5. Variation of the power of the Stirling engine for different the heat source 
temperature and x. 

solution vector X exists that minimizes all the K objective functions 
simultaneously. In other words, in various real-life problems, objec¬ 
tives at issues conflict each other. Thus, optimizingX with respect to 
a uni-objective often leads to unacceptable results with respect to 
the other objectives. Thus, a new concept, well-known as the “Par¬ 
eto optimum solution”, is used in multi-objective optimization 
problems. A feasible solution X is called “Pareto optimal” if there 
exists no other feasible solution Y that dominates solution X. By def¬ 
inition, a feasible solution Y is said to dominate another feasible 
solution X, if and only if, f(Y) ^ f(X) for i = 1,..., k and f{Y) < ffX) 
for at least one objective function. This means that a feasible vector 
X is called Pareto optimal if there is no other feasible solution Y that 
would reduce some objective function without causing a simulta¬ 
neous increase in at least one other objective function. The set of 
all feasible non-dominated solutions in X is referred to as the 
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Fig. 6. Variation of the power of the Stirling engine for different the heat source 
temperature and A. 



Fig. 7. Variation of the power of the Stirling engine for different A and UA. 


“Pareto optimal set”, and for a given Pareto optimal set, the corre¬ 
sponding objective function values in the objective space are called 
the “Pareto optimal frontier” [48 . 

For having a helpful insight into the concept of Pareto optimal 
consider Fig. 12. This figure demonstrates the objective functions 
space for an optimization problem with two objective functions 
/i and f 2 . The feasible and infeasible areas are demonstrated as 
the areas inside and outside of the curve respectively. Every point 
in the feasible area is a solution of the problem. In Fig. 12, the val¬ 
ues of both functions/i and f 2 . for point M are lower than the cor¬ 
responding values of point j. Thus, point M dominates point j or 
point M is better than point j. In the same way, points L and N dom¬ 
inate M. But points like I and I< neither dominate M nor are domi¬ 
nated by M. Therefore, only those points that are located in the left- 
down parts of M will dominate M (like L and N). Now consider 
point R , there is no point in the feasible area located in the left- 
down part of it. Thus, R is not dominated by any other feasible 
point. Therefore, R is a Pareto optimal. This is true for points 
P, A , R , E, T, 0 and all of the other points located at the bold curve 
indicated as “The Pareto Optimal Frontier” in Fig. 12. Notice that, 
in the feasible area, the minimum values of /j and f 2 belong to 



Fig. 8. Variation of the thermal efficiency of the Stirling engine for different A and x. 



Fig. 9. Variation of the maximum power of the Stirling engine for different A and x. 

points P and 0, respectively. Thus, P and 0 are the solutions for 
uni-objective optimization problems which their objective func¬ 
tions are/i and / 2 , respectively. Other points on the Pareto frontier 
are also optimal, but which of them should be chosen as the final 
solution. For doing this a decision-making process is required. 

2.2.2. Optimization via EA 

In the present study, the Pareto frontier is achieved using the 
Genetic Algorithm (GA) which is a branch of evolutionary algo¬ 
rithm. Genetic Algorithms were evolved by John Flolland in the 
1960s to import the mechanisms of natural adaptation into com¬ 
puter algorithms and numerical optimization [50]. They are used 
as a computer simulation in a way that a population of abstract 
representations (called chromosomes or the genotype of the gen¬ 
ome) of candidate solutions (called individuals, creatures, or phe- 
no-types) to an optimization problem develops toward better 
solutions. It begins from a population of randomly produced exclu¬ 
sives and occurs in generations. In every generation, the fitness of 
each individual is measured; multiple individuals are stochasti¬ 
cally chosen from the current population, and altered to form a 
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Fig. 10. Variation of the maximum power of the Stirling engine for different X and 

m)reg. 



Fig. 12. Schematic of the objectives space [35]. 


22.3. Objective functions, decision variables and constraints 

Two objective functions for this study are the system’s output 
power (P) and the Stirling engine’s thermal efficiency (^ t ) denoted 
by Eqs. (9) and (12), respectively. 

In this paper six decision variables have been considered as 
follow: 



Fig. 11. Variation of the power of the Stirling engine for different UA and x . 


new population. The new group is employed in the next iteration. 
Generally, the algorithm finishes when either a maximum number 
of generations have been produced, or a satisfactory fitness level 
has been gained for the population. If the algorithm has ended ow¬ 
ing to a maximum number of generations, a satisfactory solution 
might have been reached or might not have been accomplished. 
In genetic algorithms, a candidate solution to a problem is named 
a chromosome, and the evolutionary viability of each chromosome 
is shown by a fitness function. This method can be applied to opti¬ 
mize nonlinear problems [48,51 . 

Multi-objective evolutionary algorithms (MOEAs) have been 
evolved in the recent decades by a lot of tests on complex mathe¬ 
matical problems. On real-world it has shown that they can omit 
the obstacles of classical methods [48,51]. The structure of the 
MOEA exploited in this investigation is depicted in Fig. 13 [50]. 
The real values of decision variables are applied instead of their 
binary coded. 


T h Temperature of heat source 

A r Heat transfer area ratio 

ib Ml 

Y Oh 

T h Temperatures of the working fluid in the high 
temperature isothermal process 3-4 
x Temperature ratio 

U h A h Heat transfer coefficient between the engine and the 
high-temperature source 


The objective functions with respect to following constraints 


have been solved: 


1100 ^ T h ^ 1300 K 

(16) 

800 ^ T h ^ 1000 K 

(17) 

0.5 ^ t// ^ 1 

(18) 

0.5 ^ Ar ^ 1 

(19) 

0.5 ^ x ^ 0.7 

(20) 



Fig. 13. Scheme for the multi-objective evolutionary algorithm used in the present 
study [34-36]. 
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l .5 < i/ H a H < 6.5 i<w ic 1 


2.3. Decision-making in the multi-objective optimization 


(21) up the fuzzy sets (Eq. (26)) is turned out as a result of the intersec¬ 
tion D = nwith a membership function, 

Ho (X) = njL, H A (X) = min \i A (x) X € L (26) 

1 J= l,-,k 1 


In multi-objective optimization, a procedure of decision-mak¬ 
ing for opting the final optimal solution is implicated. There are a 
lot of decision-making methods and these methods can be utilized 
in decision-making to choose the final optimal solution from the 
Pareto frontier which is achieved by optimization. Owing to differ¬ 
ent dimensions of various objectives in a multi-objective optimiza¬ 
tion problem and necessity of unified dimensions and scales of 
objectives space, in turn, objectives vectors should be non-dimen- 
sionalized before the decision-making process. Some methods such 
as linear, Euclidian, and fuzzy can be utilized for non- 
dimensionalization. 


• Linear non-dimensionalization 


Consider the matrix of objectives at various points of the Pareto 
frontier is denoted by Fy where i is index for each point on the Par¬ 
eto frontier and j is the index for each objective in objectives space. 
Therefore, a non-dimensionalized objective, FjJ, is defined as, 


Fjj =-y— for maximizing objective 

max(Fy) 


(22a) 



-- 2 ——— for minimalizing objective 

max(l/F,j) 


(22b) 


• Euclidian non-dimensionalization 
In this method, a non-dimensionalized objective, F'J, is defined 
as, 


Fj] = . 1J . . = for minimizing and maximizing objectives 

(23) 

• Fuzzy non-dimensionalization 

In this method, a non-dimensionalized objective, Fj|, is defined 
as, 



Fjj - min(Fij) 
max(Fy) - min(Fy) 


for maximizing objectives 


(24a) 



max(Fy) - Fy 
max(Fy) - min(Fy) 


for minimizing objectives 


(24b) 


In the present research, most well-known kind of decision-making 
procedures such as the fuzzy Bellman-Zadeh, LINMAP and TOPSIS 
method are employed in parallel and the final optimal solution 
has been chosen based on real-world experience and criteria among 
the suggested solutions. The Bellman-Zadeh method uses the fuzzy 
non-dimensionalization while the other methods (LINMAP and 
TOPSIS) utilize Euclidian non-dimensionalization. The mentioned 
decision-making algorithms are represented in the following. 


2.3.1. Bellman-Zadeh decision-making method 

When using the Bellman-Zadeh approach, each F/X) is replaced 
by a fuzzy objective function or a fuzzy set, 

Aj = {X,X e L,j = 1,2,... k (25) 

where ju A .(X) is a membership function of Aj [57]. 

A final decision is defined by the Bellman and Zadeh model as 
the intersection of all fuzzy criteria and constraints and is repre¬ 
sented by its membership function. A fuzzy solution D with setting 


Using Eq. (24a), it is possible to obtain the solution proving the 
maximum degree as follows, 


max fi D (X) = max min fi A .{x) 

XgL j=1 J 


X° = arg max min p A . (x) 

XeL j=l,...,k J 


(27) 

(28) 


To obtain Eq. (27), it is necessary to build membership functions 
ju A (X) ;j = 1,... ,k, reflecting a degree of achieving “own” optima 
by the corresponding F)(X), X e L; j = 1,..., k. This is satisfied by the 
use of the membership functions [58]. The membership function 
of objectives and constraints, linear or nonlinear, can be chosen 
depending on the context of problem. One of possible fuzzy convo¬ 
lution schemes is presented below [59 . 

Initial approximation for X-vector is chosen. Maximum (mini¬ 
mum) values for each criterion F/X) are established via scalar max¬ 
imization (minimization). Results are denoted as “ideal” points 
{Xj°J= 

The matrix table {F}, where the diagonal elements are “ideal” 
points, is defined as follows: 



F ,(X?) F 2 (X°)... F„(X?) 
F,(X°) F 2 (X°)... F„(X°) 


L F i (X^) F 2 (X° n )... F„(X°)J 


(29) 


Maximum and minimum bounds for the criteria are defined: 


F™ n = min F,(X°) i=l,...,n 
j J 

(30a) 

F| mx = max Fj(X°) i= 
j 

(30b) 


The membership functions are assumed for all fuzzy goals as fol¬ 
lows.For minimized objective functions: 


Mr(X) = 


o if f,(x) > f; 

pmax_p 


max 


pmax _pmin 

1 


if F™ n <F, ^Fr x 
if Fi(x)^Ff 


I N 1 ! 
:min 


For maximized objective functions: 


Mh(X) = 


p._pmin 
pmax _pm\n 

0 


-max 


if Fi(x) > F" 
if Ff n < F, 5 F™ x 

if F,-(x) < F" 


I N 1 ! 
-min 


Fuzzy constraints are formulated: 

Cj(X) < G J max + dj, j = 1,2,...,k 


(31a) 


(31b) 


(32) 


where dj is a subjective parameter that denotes a distance of 
admissible displacement for the bound G™ ax of the jth constraint? 
Corresponding membership functions are defined in following 
manner: 



Mq(X) = < 1 



if Gfx) > G™ x 

/-max 

if GJ"“ < G,(X) < G, max + dj 
if Gfx) < G™ x 
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thermal efficiency 


Fig. 14. Pareto optimal frontier in objectives’ space. 


ideal point is selected as a final desired optimal solution, hence, i in¬ 
dex for a final solution, i fina i is, 

ifi„ a i = i e min(d i+ ); i=l,2,...,m (36) 


2.3.3. TOPSIS decision-making method 

In this method, beside the ideal point a non-ideal point is de¬ 
fined. The non-ideal point is the ordinate in objectives space in 
which each objective has its worst value. Therefore, beside the 
solution distance from ideal point, d I+ , the solution distance from 
non-ideal point denoted by d,_ is used as a criterion for selection 
of the final solution. Hence, 





r;non-ideal\ 

F J ) 


2 



In continuing the TOPSIS method a parameter is defined as follows, 



d/_ 

di+ + d[_ 



In TOPSIS method, a solution with minimum C7, is selected as a de¬ 
sired final solution, therefore, if ifi na i is index for the final selected 
solution, we have, 


A final decision is determined as the intersection of all fuzzy criteria 
and constraints represented by its membership functions. This 
problem is reduced to the standard nonlinear programming prob¬ 
lems: to find such values of X and k that maximize k subject to 


ifinai = i G maxfC/i); i = l,2,...,m (39) 

2.4. Second scenario results 


A < li F i = 1,2,.. .,n 

(34) 

2^H Cj , j = 1,2,...,k 

The solution of the multi-criteria problem discloses the meaning of 
the optimality operator and depends on the decision-makers expe¬ 
rience and problem understanding. 

2.3.2. LINMAP decision-making method 

An ideal point on the Pareto frontier is the point in which each 
objective is optimized regardless to satisfaction of other objectives. 
It is clear that in the multi-objective optimization it is impossible 
to have each objective in its optimal condition that it can obtain 
in a single-objective optimization. Therefore, the ideal point is 
not located on the Pareto frontier. In the LINMAP method, after 
Euclidian non-dimensionalization of all objectives, the distance of 
each solution on the Pareto frontier from the ideal point denoted 
by d,+ is determined as follow, 

4+ = /ElL/ k-C' (35) 

where denotes the number of objective while stand for each solu¬ 
tion on the Pareto frontier (i = 1, 2,..., m). In Eq. (35), Ff eal is the 
ideal value for jth objective obtained in a single-objective optimiza¬ 
tion. In LINMAP method, the solution with minimum distance from 


The system’s output power and thermal efficiency of the Stirling 
engine are maximized simultaneously using the multi-objective 
optimizing method which works based on the NSGA-II algorithm. 

In this regard, optimization is performed with objective func¬ 
tions that are expressed by Eqs. (9) and (12) constraints which 
are expressed with Eqs. (16)-(21). Design parameters (decision 
variables) of optimization are the temperature of heat source, heat 
transfer area ratio, i/q temperatures of the working fluid in the high 
temperature isothermal process 3-4, temperature ratio, Heat 
transfer coefficient between the engine and the high-temperature 
source. 

Fig. 14 illustrates the Pareto frontier in the proposed objectives’ 
space obtained using the multi-objective evolutionary algorithm 
(NSGA-II). Three final solutions selected by the Bellman-Zadeh, 
LINMAP and TOPSIS decision-makers are indicated in Fig. 14. 

The results of two-objective optimization based on Fuzzy, LIN¬ 
MAP and TOPSIS decision making methods are shown in Table 1. 
In this Table, results for P and t] t objective functions based on six 
decision variables are categorized. Temperature of heat source 
(T h ) is obtained between 1100 I< which is inacceptable value due 
to previous studies [42]. Temperature ratio (x) is obtained between 
0.500 and 0.501 which is inacceptable range due to previous stud¬ 
ies [42]. 

Table 2 compares optimal results obtained in this paper by two 
scenarios. 


Table 1 

Decision making of multi-objective optimal solutions. 


Decision making method 

Decision variables 





Objective functions 

Th (K) 

r„(K) 


Ar 

X 

ua (kwr 1 ) 

P (kW) 


Bellman-Zadeh 

1100.7 

999.75 

0.996 

0.996 

0.501 

4.427 

77.114 

17.25 

LINMAP 

1100.6 

999.8 

0.998 

0.997 

0.501 

3.485 

68.874 

19.6 

TOPSIS 

1100.6 

999.8 

0.998 

0.996 

0.500 

1.961 

49.510 

25.04 
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Table 2 

Values of variables and objective functions obtained in two optimization scenarios. 
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Decision making method 

Variables 






Objective functions 

Th (K) 

T„ 00 


Ar 

X 

ua (kwr 1 ) 

P (kW) 


Scenario 1 

— 

1300 

970 

1.000 

1.000 

0.500 

2.000 

75.306 

17.82 

Scenario 2 

Bellman-Zadeh 

1100.7 

999.75 

0.996 

0.996 

0.501 

4.427 

77.114 

17.25 


LINMAP 

1100.6 

999.8 

0.998 

0.997 

0.501 

3.485 

68.874 

19.6 


TOPSIS 

1100.6 

999.8 

0.998 

0.996 

0.500 

1.961 

49.510 

25.04 


3. Conclusions 

In the presented paper, the effects of some parameters such as 
UA, ( UA) reg , heat source temperature (T H ) and volumetric ratio (2) 
on the output power and the thermal efficiency is investigated. It 
can be concluded that the effects of UA is more than ( UA) reg on 
the output power. The temperature of heat source is one of the 
important parameters affecting the output power and by increas¬ 
ing of temperature the power increases too. This behavior varies 
with changes in other parameters. 

The volumetric ratio has an effective role in output power and 
by increasing it, the power increases, but the slope of the curve de¬ 
creases by increasing temperature ratio; Also, the effects of 2 in the 
range 1.5-3.5 is sensible. Volumetric ratio affects thermal effi¬ 
ciency similar to the output power. Also, the impacts of x on the 
power and thermal efficiency are evaluated. The results of two sce¬ 
narios show that the theory can be useful for design and the perfor¬ 
mance analysis of Stirling heat engine. 
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